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Abstract. It is proven that the sum rules for X-ray magnetic dichroism (XMCD) spectra that are used 
to separate spin and orbital contributions to the magnetic moment are formally correct for an arbitrary 
strength of electron-electron interactions. However, their practical application for strongly correlated sys- 
tems can become complicated due to the spectral density weight spreading over a broad energy interval. 
Relevance of incoherent spectral density for the XMCD sum rules is illustrated by a simple model of a 
ferromagnet with orbital degrees of freedom. 

PACS. 78. 70. Dm X-ray absorption spectra - 78. 70. En X-ray emission spectra and fluorescence - 75. 30. Mb 
Valence fluctuation, Kondo lattice, and heavy-fermion phenomena - 71.28.+d Narrow-band systems; 
intermediate- valence solids 



X-ray magnetic circular dichroism (XMCD) [T|P] is a 
powerful technique to investigate both bulk and surface 
magnetic properties of materials. In particular, it allows 
to measure separately spin and orbital contributions to 
the magnetic moment of ferromagnets. Examples of nu- 
merous applications of this method are recent studies of 
magnetism in thin metallic films [3], cobalt nanoparticles 
and clusters 0|, in magnetite Fe3C>4 and in dilute mag- 
netic alloys [0] ■ A concrete way to separate the orbital and 
spin magnetic moments is using the XMCD sum rules [3 
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where Nh is the number of holes in d-band, /il 2 3 ( w ) are 
spectral intensities for 1^,3 spectra, A/i is the difference 
between the spectra for left and right circularly polar- 
ized radiation and /j, tot is the total absorption intensity 
for unpolarized one, S z , L z and T z are projections of the 
total spin, orbital moment and the spin dipole operator 
on the magnetization direction: S = ^ s iiI J = ^li,T = 

i i 

^2 ( s i — 3i"i (rjSj) frf) with and 1^ being the spin and 

i 

orbital moments for i-th electron, the coordinate oper- 
ators. 

A simple derivation of the XMCD sum rules in the in- 
dependent electron approximation was presented in Ref. [5j • 
This derivation is based on a purely band picture of elec- 



tron states in solids. The opposite case of strongly local- 
ized electrons which are characterized by atomic states 
with well-defined term and multiplet structure has been 
considered in Refs. [7lll()| (a somewhat more simple deriva- 
tion with the use of second quantization formalism for 
the atomic states was presented in the paper [S]). At the 
same time, many interesting systems such as magnetite 
[5] should be considered as strongly correlated systems 
demonstrating simultaneously both itinerant and local- 
ized features of "magnetic" electrons | 1 11112) . Actually, 
the sum rules are widely exploited by experimentalists 
for such systems as well. In this work we present a for- 
mal justification of this simple way to proceed. At the 
same time, we discuss separately contributions to the sum 
rules from coherent and incoherent parts of the electron 
spectral density. We demonstrate that a proper account 
of the incoherent (nonquasiparticle) contributions is nec- 
essary for consistent treatment of the XMCD spectra of 
strongly correlated systems. 

The spectral intensity (J,(u>) of the X-ray absorption 
and emission spectra (XAS and XES, respectively) is de- 
termined in the dipole approximation by the imaginary 
part of the Green's function lj 

((p-e*|p.e)) ti; = e*e^H (2) 

OO 

G a p M = ((P.ML = -% f dte^ (\p a (t) , Pfi \) (3) 

o 

where e is the photon polarization vector and p is the 
momentum operator, a, [3 are the Cartesian indices and 
the brackets stand for the Gibbs average in the initial 
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state. For the case of XAS only the part of this opera- 
tor works which corresponds to the transitions from core 
states | a) (with the annihilation operators b a ) to the con- 
duction electron states |A) created by operator c\ 

pi +) =J2(Mp a \a)c{b a: (4) 

a A 

and only the Hermitian conjugated part pit works for 
the XES. For the case of 1,2(1^3) spectra a labels total mo- 
ment projection for 2p l / 2 (2p 3 ^ 2 ) states, correspondingly, 
and A labels spin projection and orbital indices for 3d- 
electrons. Transitions to s-states, which are also allowed in 
the dipole approximation, are irrelevant for magnetism (in 
particular, they practically do not contribute to XMCD) 
and therefore will be neglected further. 

It is important that Eqs.|0, @ are formally exact ir- 
respective to the degree of localization or derealization 
of d-electrons. The use of the atomic representation for 
the states |A) allows to obtain explicit expressions for the 
matrix elements (X\p a \&) in terms of 3?ij-symbols, frac- 
tional parentage coefficients and irreducible matrix ele- 
ments 0IE1E3]- However, for the derivation of the sum 
rules such a concretization is not necessary and using the 
A-representation yields a simpler way to proceed. Actually, 
total intensities of the XAS for different photon polariza- 
tions which are needed to obtain the XMCD sum rules are 
determined by the integrals of the spectral density in the 
infinite limits. Due to the Kramers-Kronig relations, the 
latter are equal to 

oo 

/d 
— lmGa/3 M = hm [uiReG a p (u)] = ([p a ,pp]) 
7T u— »oo 

— oo 

(5) 

Here we have used the equations of motion for the Green's 
functions On substituting Eq.(0J into the right-hand 
side of Eq.JSJ) one can see that the commutator contains 
just the one-particle density matrices which are expressed 
in terms of the corresponding anticommutator Green's 
functions: 

oo 

Pva = (cW) = - J ^f(E)lm((cy\4)) E (6) 

— oo 

where f(E) is the Fermi function; similarly we introduce 
the core density matrix 

(%? = Hb a r) (7) 

Assuming that in the initial state the core electron states 
are completely occupied, i.e., p™ r , e = S aa /, one can receive 
the expression for the total spectral weight (integrated 
intensity) of h 2 ^ spectra: 

M? = 27r (i m il' m ^°')^ m 2 <7 IP' e * ( 8 ) 

rrijinarn' a' 

x(l- p)p ■ e\lm'-a'){lm'-a'\jmj) 



which differs from Eq.(2) of Ref. [H] only by the replace- 
ment of one-electron expression for the density matrix 
pmm' by the exact one. Here I = 2, j = 1/2 (3/2) for 
L2 (L3) spectra, m,m' and a, a' are the orbital and spin 
projections of d-electrons, respectively, rrij is the total 
moment projection for the core states. It is worthwhile 
to stress that using the orbital indices for itinerant elec- 
trons does not mean any approximation: for any partic- 
ular method of band structure calculations it is always 
possible to re-expand the Wannier functions at a given 
site into the spherical harmonics. Thus the density ma- 
trix pmm' is in general a linear combination of the band 
structure occupation numbers with a proper symmetry. 
Since the total spin and orbital moment as well as the spin 
dipole operator are one-particle operators, their averages 
are completely determined by the density matrix. Further 
use of the Wigner-Eckart theorem to extract angular de- 
pendences of the matrix elements and transformations of 
the arising 3j-symbol products repeat the derivation in 
Ref. 9 . Therefore the XMCD sum rules are formally 
valid without any restrictions. 

On the other hand, the values of (S z ) and (L z ) ob- 
tained in band calculations can violate the sum rules. The 
calculation of the one-particle density matrix for strongly 
correlated systems remains a quite nontrivial problem. 
Generally, this quantity contains both coherent (quasipar- 
ticle) contributions which are formally connected with the 
poles of the Green's function and incoherent (nonquasipar- 
ticle) ones which are formally connected with the branch 
cuts |15j . In a number of cases (for example, for strongly 
correlated metals in the vicinity of the Mott transition or 
for doped Mott insulators ^B]) the quasiparticle spectral 
density is much narrower than the incoherent one. One 
should be therefore careful to avoid a confusion of the in- 
coherent part with a spectrum background. 

In practice, the use of the XMCD sum rules for strongly 
correlated systems is not a completely well-defined proce- 
dure since it is impossible to integrate the spectral density 
in infinite limits. Atomic multiplet structure leads some- 
times to much broader distribution of the spectral den- 
sity in comparison with a standard band picture. This 
fact has been recently demonstrated ^2] for the case of 
(La,Sr)Mn03 where the configuration-interaction calcu- 
lations for Mn 3+ ion gave the spectral density with the 
width of order of 6-^7 eV, in comparison with the value 
2-^3 eV within the local density approximation (LDA) or 
LDA+U. At the same time, the energy separation between 
L2 and L3 spectra varies from 6-^8 eV for the light 3d met- 
als (Ti, V) to 15^20 eV for heavy ones (Co, Ni, and Cu). 
This means that for strongly correlated systems L2 and 
L 3 spectra can overlap appreciably. The energy distribu- 
tion of the atomic multiplet structure can be comparable 
with the spin-orbit splitting of the relevant core levels also 
for rare earth systems which leads to some problems with 
the practical applications of the XMCD sum rules in the 
latter case [T%] . 

It is a common practice to interpret the core-level spec- 
tra for isolated atoms in terms of many-electron multi- 
plet picture rather than one-electron quantum numbers 
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| A) used above. These two approaches can be related via 
the representation of one-electron operators in terms of 
the Hubbard X-operators QUI X r < r ' = \r)(r'\ where P = 
{nLSMU} labels atomic configuration n, term LS and 
moment projections M, E: 

4ma= £ s/nGrl^LnM^Ln-iMn-xlm) 

x(S n E n \S n -xE n ^a)X r -< r ^ (9) 

where Gp" _ are the fractional parentage coefficients (see 
Refs. [20121) V Eq.© enables one to reproduce the sum 
rules in the form yielding a detailed information concern- 
ing term and multiplet structure [7II511TU] . 

At the same time, for a solid the atomic description 
can be inappropriate. In particular, ferromagnetism itself 
is an essentially band phenomenon. Moreover, it cannot be 
properly described in the simplest Hubbard-I approxima- 
tion [22j which assumes a formation of individual Hubbard 
subbands from separate transitions between the atomic 
levels. In particular, for a narrow-band ferromagnet we 
cannot satisfy the sum rules (kinematic relations) for the 
X-operators in this approximation. 

To demonstrate this we consider a simple model of 
a narrow-band itinerant electron ferromagnet with orbital 
degrees of freedom which generalizes the standard infinite- 
U Hubbard model: 

n = t*r n x°jrxr , ° (io) 

kmcr 

where t km is the orbital-dependent band energy, X^'^ 3 is 
the Fourier transform of the Hubbard operators and 
labels the hole state at a site. 

For this model the exact sum rules should be satisfied 
for arbitrary to, a: 

r /V'0,0\ / yO,OTn yinn.Oi \ 1 / v am v am,0\ 

o = n = {X ) = {X i \ l ) = 2^\ A -k A k / 

k 

+ 00 

= -£ / —f(E)M(xr'°\x°jr))E. (11) 

K -00 

where S is the concentration of current carriers (holes). 
The Hubbard-I approximation for this model reads 

((X° m -°\xX m ))E = [E- t km (N + N am )}-\ (12) 

According to this expression, the quasiparticle pole for 
a =\ corresponds to a narrowed band and lies (for the 
saturated ferromagnet case) above the Fermi level of the 
holes which obviously violates the sum rule 1)1 ]|. In fact, 
these sum rules are satisfied only due to incoherent (non- 
quasiparticle) states which are present below the Fermi 
level for the Hubbard model Similar to these papers 
one obtains for the minority-spin Green's functions in the 
leading order in small parameter S 

- hm((Xi m ' \X°^ m )) E 



). (13) 

qm' 

where E^ m = ikm(-Wo +N-[ m ), w qm m' are the frequencies 
of the corresponding spin-flip transitions. Thus we have 
100% incoherent contribution to the spectral density be- 
low the Fermi level. As the hole concentration increases, 
the ferromagnetic state becomes non-saturated and a nar- 
row quasiparticle minority-spin band occurs |24j . How- 
ever, the main spectral density for this spin projection is 
still due to the nonquasiparticle contribution. This exam- 
ple demonstrates that, despite the XMCD sum rules for 
strongly correlated ferromagnets have formally a standard 
"one-electron" form, the energy distribution of the spec- 
tral weight can be drastically different from usual band 
picture. 

The theoretical consideration shows that the non- 
quasiparticle contributions to XES and XAS should be 
clearly observed. Moreover, such contributions should be 
considerably enhanced by the interaction with the core 
hole As well as for XES and XAS in general, the 
main nonquasiparticle effects in the XMCD are connected 
with the occurrence of the incoherent spectral density in 
the energy gap for the half-metallic ferromagnets |2b|. 
These states arise either only below the Fermi energy 
(for the majority-gap half-metallic ferromagnets such as 
magnetite) or only above it (for the minority-gap half- 
metallic ferromagnets such as Heusler alloys or CrO^); 
therefore they can be studied by XES and XAS meth- 
ods, correspondingly. Whereas a standard band theory 
predicts 100% polarization for the conduction electron (or 
hole) states, the depolarization due to the nonquasiparti- 
cle states can be very strong (for example, in the infinite-?/ 
Hubbard model limit there is no polarization at all [23)). 
At the same time, ab initio calculations of the correlation 
effects for the half-metallic Heusler alloy NiMnSb give a 
rather small spectral weight of the nonquasiparticle states 
(about 4%) [27] and therefore it would be preferable to 
investigate them for the half-metallic ferromagnets with 
more strong correlations such as Fe3C>4 (by XES) and 
CrOa (by XAS). 
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